f{R) gravity constrained by PPN parameters and stochastic background of 

gravitational waves 



S. Capozziello\ M. De Laurentis2\ S. Nojiri^, S. D. Odintsov"* 
^ Dipariimento di Scienze fisiche, Universita di Napoli " Federico II", INFN Sez. di Napoli, 
Compl. Univ. di Monte S. Angela, Edificio G, Via Cinthia, 1-80126, Napoli, Italy 
^ Politecnico di Torino and INFN Sez. di Torino, 
Corso Duca degli Ahruzzi 24, 1-10129 Torino, Italy 
^ Department of Physics, Nagoya University, Nagoya 4S4-8602, Japan 
'^Institucio Catalana de Recerca i Estudis Avancats (ICREA) and Institut de Ciencies de I Espai (lEEC-CSIC) , 
Campus UAB, Facultat de Ciencies, Torre C5-Par-2a pi, E-08193 Bellaterra (Barcelona), Spain. 

(Dated: August 11, 2008) 

We analyze seven different viable /(7?)-gravities towards the Solar System tests and stochastic 
gravitational waves background. The aim is to achieve experimental bounds for the theory at local 
and cosmological scales in order to select models capable of addressing the accelerating cosmological 
expansion without cosmological constant but evading the weak field constraints. Beside large scale 
structure and galactic dynamics, these bounds can be considered complimentary in order to select 
self-consistent theories of gravity working at the infrared limit. It is demonstrated that seven viable 
/(ii)-gravities under consideration not only satisfy the local tests, but additionally, pass the above 
PPN-and stochastic gravitational waves bounds for large classes of parameters. 
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I. INTRODUCTION 



The currently observed accelerated expansion of the Universe suggests that cosmic flow dynamics is dominated by 
some unknown form of dark energy characterized by a large negative pressure. This picture comes out when such 
a new ingredient, beside baryonic and dark matter, is considered as a source in the r.h.s. of the fleld equations. 
Essentially, it should be some form of un-clustered, non-zero vacuum energy which, together with (clustered) dark 
matter, should drive the global cosmic dynamics. 

Among the proposals to explain the experimental situation, the "concordance model", addressed as ACDM, gives a 
reliable snapshot of the today observed Universe according to the CMBR, LSS and SNela data, but presents dramatic 
shortcomings as the ''coincidence and cosmological constant problems" which point out its inadequacy to fully trace 
back the cosmological dynamics jli]. 

On the other hand, alternative theories of gravity, extending in some way General Relativity (GR), allows to 
pursue a different approach giving rise to suitable cosmological models where a late-time accelerated expansion can 
be achieved in several ways. This viewpoint does not require to flnd out candidates for dark energy and dark matter 
at fundamental level (they have not been detected up to now), it takes into account only the "observed" ingredients 
(i.e. gravity, radiation and baryonic matter), but the l.h.s. of the Einstein equations has to be modified. Despite of 
this modification, it could be in agreement with the spirit of GR since the only request is that the Hilbert-Einstein 
action should be generalized asking for a gravitational interaction acting, in principle, in different ways at different 
scales 

The idea that Einstein gravity should be extended or corrected at large scales (infrared limit) or at high ener- 
gies (ultraviolet limit) is suggested by several theoretical and observational issues. Quantum field theory in curved 
spacetimes, as well as the low-energy limit of String/M theory, both imply semi-classical effective actions containing 
higher-order curvature invariants or scalar-tensor terms. In addition, GR has been definitely tested only at Solar 
System scales while it may show several shortcomings if checked at higher energies or larger scales. Besides, the 
Solar System experiments are, up to now, not so conclusive to state that the only viable theory ofgravity is GR: for 
example, the Hmits on PPN parameters should be greatly improved to fully remove degeneracies 

Of course, modifying the gravitational action asks for several fundamental challenges. These models can exhibit 
instabilities 0| or ghost -like behavior while, on the other hand, they have to be matched with observations and 
experiments in the appropriate low energy limit. 

Despite of all these issues, in the last years, some interesting results have been achieved in the framework of the so 
called /(i?)-gravity at cosmological. Galactic and Solar System scales. Here f{R) is a general (analytic) function of 
the Ricci scalar R (see Refs. [Q, l3, iS] for review). 
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For example, there exist cosmological solutions that give the accelerated expansion of the universe at late times 
0, [13, [m, [l2|- In addition, it has been discovered that some stability conditions can lead to avoid ghost and tachyon 
solutions. Furthermore there exist viable .f(R) models which satisfy both background cosmological constraints and 
stability conditions [TtI. fisl. fisl. 1 2d [21I and results have been achieved in order to place constraints on f{R) 

cosmological models by CMBR anisotropies and galaxy power spectrum (2^. [23l . [23l . Moreover, some of such viable 
models lead to the unification of early-time inflation with late-time acceleration [l9l. [2ol. |2]]|. 

On the other hand, by considering /(i?)-gravity in the low energy Hmit, it is possible to obtain corrected gravitational 
potentials capable of explaining the flat rotation curves of spiral galaxies or the dynamics of galaxy clusters without 
considering huge amounts of dark matter [H, [2^, ^ [H, [li, [s^l • 

Furthermore^several authors have dealt with the weak field limit of fourth order gravity, in particular considering 
the PPN limit [13, [13, H [H, [H, [13, [H, Sil and the spherically symmetric solutions [13, [3, [U, [ii] . 

This great deal of work needs an essential issue to be pursued: we need to compare experiments and probes at local 
scales (e.g. Solar System) with experiments and probes at large scales (Galaxy, extragalactic scales, cosmology) in 
order to achieve self-consistent f{R) models. Some work has been done in this direction (see e.g. [la|) but the large 
part of efforts has been devoted to address single data sets (observations at a given redshift) by a single model which, 
several time, is not working at other scales than the one considered. In particular, a given f{R) model, evading Solar 
System tests, should be not simply extrapolated at extragalactic and cosmological scales only requiring accelerated 
cosmological solutions but it should be confronted with data and probes coming from cosmological observations. 
Reliable models are then those matching data at very different scales (and redshifts). 

In order to constrain further viable /(-/?)-models, one could take into account also the stochastic background of 
gravitational waves (GW) which, together with cosmic microwave background radiation (CMBR), would carry a 
huge amount of information on the early stages of the Universe evolution. In fact, if detected, such a background 
could constitute a further probe for these theories at very high red-shift [g^I- On the other hand, a key role for the 
production and the detection of the relic gravitational radiation background is played by the adopted theory of gravity 
[H, [i3|- This means that the effective theory of gravity should be probed at zero, intermediate and high redshifts to 
be consistent at all scales and not simply extrapolated up to the last scattering surface, as in the case of GR. 

The aim of this paper is to discuss the PPN Solar-System constraints and the G W stochastic background considering 
some recently proposed f{R) gravity models (l3l . \l4 . [13, [l3, [13, [HI which satisfy both cosmological and stability 
conditions mentioned above. Using the definition of PPN-parameters 7 and (3 in terms of /(i?)-models and the 
definition of scalar GWs [i^l, we compare and discuss if it is possible to search for parameter ranges of /(i?)-models 
working at Solar System and GW stochastic background scale. This phenomenological approach is complementary to 
the one proposed, e.g. in [13, [l^l where also galactic and cosmological scales have been considered to constraint the 
models. 

The layout of the paper is the following. In Sec. II, we review the field equations of f{R) gravity in the metric 
approach and their scalar-tensor representation, useful to compare the theory with observations. In Sec. Ill, we review 
and discuss some viable f{R) models capable of satisfying both local gravity prescriptions as well as the observed 
cosmological behavior. In particular, we discuss their stability conditions and the field values which have to achieved 
to fulfill physical bounds. Sec. IV is devoted to derive the values of model parameters in agreement with the PPN 
experimental constraints while, in Sec. V, we deal with the constraints coming from the stochastic background of GWs. 
These latter ones have to be confronted with those coming from PPN parameterization. Discussion and conclusions 
are drawn in Sec. VI. As a general remark, we find out that bounds coming from the interferometric ground-based 
(VIRGO, LIGO) and space (LISA) experiments could constitute a further probe for f{R) gravity if matched with 
bounds at other scales. 



II. f{R) GRAVITY 

Let us start from the following action 

S ^Sg+Sm = ^ j d^x^[R + f{R) + , (1) 

where we have considered the gravitational and matter contributions and k"^ = IGttG. The non-linear f{R) term 
has been put in evidence with respect to the standard Hilbert-Einstein term R and Cm is the perfect-fluid matter 
Lagrangian. The field equations are 

i9^.i^(i?) - R^.F'iR) - g^.DF'iR) + V^V.F'(i?) - -^Tj^Z'l (2) 
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Here F{R) = R + f{R) and T/;™^s the matter energy - momentum tensor. By introducing the auxiliary field A, one 
can rewrite the gravitational part in the Action ^ as 



Sa = ^j d'xV^ {(1 + fiA)) {R-A)+A + /(A)} . (3) 

As it is clear from Eq.(l3]), if F'{R) = 1 + f'{R) < 0, the coupling k^jj — k"^ /F'{A) becomes negative and the theory 
enters the anti-gravity regime. Note that it is not the case for the standard GR. 

Action ^ can be recast in a scalar-tensor form. By using the conformal scale transformation g^^ e'^gtj.L' with 
(T = — In (1 + f'{A)), the action can be written in the Einstein frame as follows 0|: 



Se = 



'dpOdaO - y(cr) 



(4) 



where 



V[o) = e^g (e--) - e^/ (, (e-^)) = ^ - (5) 

The form of gifi""^) is given by solving a = — In (1 -t- /'(A)) = InF'(^) as ^ = g{e^°'). The transformation 
9(11^ e^9fiv induces a coupling of the scalar field a with matter. 
In general, an effective mass for a is defined as [2l| 



m„ 



1 (PV{(j) 

2 da^ 



A 



AF{A) 



1 



F'{A) {F'{A)y F"{A) 



(6) 



which, in the weak field limit, could induce corrections to the Newton law. This allows, as it is well known, to deal with 
the extra degrees of freedom of /(i?)-gravity as an effective scalar field which reveals particularly useful in considering 
"chameleon" models This "parameterization" will be particularly useful to deal with the scalar component of 

GWs. 



III. f{R) VIABLE MODELS 

Let us consider now a class of f{R) models which do not contain cosmological constant and are explicitly designed to 
satisfy cosmological and Solar-System constraints in given limits of the parameter space. In practice, we choose a class 
of functional forms of f{R) capable of matching, in principle, observational data (see for the general approach). 
Firstly, the cosmological model should reproduce the CMBR constraints in the high-redshift regime (which agree with 
the presence of an effective cosmological constant). Secondly, it should give rise to an accelerated expansion, at low 
redshift, according to the ACDM model. Thirdly, there should be sufficient degrees of freedom in theparameterization 
to encompass low redshift phenomena (e.g. the large scale structure) according to the observations [2J|. Finally, small 
deviations from GR should be consistent with Solar System tests. All these requirements suggest that we can assume 
the limits 

lim f[R) — constant, (7) 



^im/(i?) = 0, (8) 
which are satisfied by a general class of broken power law models, proposed in which are 



fi{R) = -rr? 



C2 



(9) 



or otherwise written as 
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ny^) = ^ixc^^ (10) 



where m is a mass scale and Ci_2 are dimensionless parameters. 
Besides, another viable class of models was proposed in 



Fu{R) = R + \Rc 



(11) 



Since F{R = 0) = 0, the cosmological constant has to disappear in a flat spacetime. The parameters {n, p, A, Rc} 
are constants which should be determined by experimental bounds. 

Other interesting models with similar features have been studied in (itI . [T§ . [l9l . [20l . [2l| . In all these models, a 
de-Sitter stability point, responsible for the late-time acceleration, exists for R = Ri (> 0), where Ri is derived by 
solving the equation i?i/__R(i?i) = 2/(_Ri) [i^l- For example, in the model (fTTj) . we have Ri/Rc = 3.38 for A = 2 and 
p = 1. If A is of the unit order, Ri is of the same order of Rc. The stability conditions, f^R > and /^rr > 0, are 
fulfilled for R > Ri (l3l. Moreover the models satisfy the conditions for the cosmological viability that gives rise 
to the sequence of radiation, matter and accelerated epochs 

In the region R^ Rc both classes of models ^ and (flTI) behave as 



Fiii{R)~R~\Rc 1-iRc/R) 



•,2s 



(12) 



where s is a positive constant. The model approaches ACDM in the limit R/ Rc 
Finally, let also consider the class of models [3, [H, HH 



oo. 



Fiv{R) = R-\Rc(^^y . (13) 

Also in this case A, q and Rc are positive constants (note that n, p, s and q have to converge toward the same values to 
match the observations). We do not consider the models whit negative q, because they suffer for instability problems 

associated with negative Frr [13, S|- In Fig.|(l|), we have plotted some of the selected models as function of — for 

Rc 

suitable values of {p, n, q, s. A} . 

Let us now estimate nia- for the models discussed above. For Model I [l6|, when the curvature is large, we find 

+ ^ + (14) 



and obtain the following expression: 



9 9 / n \ n + 2 



" 2n{n + l)ci ' ^ ' 



Here the order of the mass-dimensional parameter rr? should be rr? ^ 10 ^^eV^. Then in Solar System, where 
R - lO-f^^eV^ the mass is given by ~ 10 ss-i-Sn .^^j^jjg ^j^g Earth atmosphere, where i? 10 ^"eV^, it 

has to be ml ~ 10"^^+^'*" eV^. The order of the radius of the Earth is 10^ m ~ (lO'^'^eV) ^. Therefore the scalar 
field a is enough heavy if n 1 and the correction to the Newton law is not observed, being extremely small. In fact, 
if we choose n = 10, the order of the Compton length of the scalar field a becomes that of the Earth radius. On the 
other hand, in the Earth atmosphere, if we choose n = 10, for example, we find that the mass is extremely large: 

^ 10^3 Ggv _ 1029 X Mpianck ■ (16) 



Here Mpianck is the Planck mass. Hence, the Newton law correction should be extremely small. 
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- ■ - ■ I model 

II model 

Ill model 

IV model 

f(R)=R 




Figure 1: Plots of four different F{R) models as function of Model I in Eq. ^ with n = 1 and A = 2 (dashed line). Model 
II in Eq.(ITT|) with p = 2, A = 0.95 (dashdot line). Model III in Eq.(Il2l) with s = 0.5 and A = 1.5 (dotted). Model IV in Ea. lfTSll 
with q — 0.5 and A = 0.5 (solid line). We also plot F{R) = R (solid thick line) to see whether or not the stability condition 
F,R > is violated. 



In Model II 



fii{R) 



-XRo 



(17) 



if R is large compared with Rq, whose order of magnitude is that of the curvature in the present universe, we find 



fii{R) 



'XRq + X- 



R^^ 
R^P 



+ 



(18) 



By comparing Eq. (fT8|) with Eq. lfT4)) . if the curvature is large enough when compared with Rq or m^, as in the Solar 
System or on the Earth, we can set the following identifications: 



XRf) 



m?ci 



XRl^ 



+1 



2p 4-> n 



(19) 



C2 C2 

We have 41to^ ~ i?o- Then, if p is large enough, there is no correction to the Newton law as in Model I given by 
Eq.lfTO|. 

Let us now discuss the instability of fluid matter proposed in [i^, which may appear if the matter-energy density 
(or the scalar curvature) is large enough when compared with the average density the Universe, as it is inside the 
Earth. Considering the trace of the above field equations and with a little algebra, one obtains 

.2 



F^'HR) ^ , F'{R)R 2F{R) 

F(2)(i?) P 3i^(2)(i?) 3F(2)(i?) 6F(^\R)^ 
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Here T is the trace of the matter energy- momentum tensor: T = Tp'"-*''. We also denote the derivative d"'F{R)/dR'^ 
by F^"\R). Let us now consider the perturbation of the Einstein gravity solutions. We denote the scalar curvature, 
given by the matter density in the Einstein gravity, by Rb ~ {n^ /2)p > and separate the scalar curvature R into 
the sum of Rb (background) and the perturbed part Rp as R = Rb + Rp {\Rp\ ^ \Rb\)- Then Eg . (1201) leads to the 
perturbed equation: 

n - nfi- +£^!lMv RVPR F'iRb)Rb 
° - ^^'^Fi^)^'^'^ ^'^ 3F(^){Rb) 

^^^^'^ ' DRp + 2^^^WpRbW''Rp + U{Rb)Rp . (21) 



3F(^){Rb) 3F(^){Rb) " F(2)(i?, 

Here the potential U{Rb) is given by 



FW(Rb)F^^HRb)Rb F(i)(i?b) , 2F(i?b)F(3)(i?,) F(^HRb)Rb 



3^^(2)(i?fc)2 3F(2)(i?fc) 3F(^){Rb)^ 3Fi^){Rb)^ 



(22) 



It is convenient to consider the case where Rb and Rp are uniform and do not depend on the spatial coordinates. Hence, 
the d'Alembert operator can be replaced by the second derivative with respect to the time, that is: ClRp —dfRp. 
Eq. l|22p assumes the following structure: 

^ -d^Rp + U{Rb)Rp + const . (23) 



If U{Rb) > 0, Rp becomes exponentially large with time, i.e. Rp ~ ev^^^*, and the system becomes unstable. 
In the l/i?-model, considering the background values, we find 



U{Rb) = -Rb + -\^ ^ ^ {10 



sec 



G/i"* \gcm 
Rb^{lOhec)-'(^^) . (24) 



g cm ^ 



Here the mass parameter /i is of the order 



lO^^sec - (lO^^^eV) ^ . (25) 



Eq. l|24p tells us that the model is unstable and it would decay in 10 ^6 ggg (considering the Earth size). In Model I, 
however, U{Rb) is negative: 

U{R^) ~ r^— f < . 26 

cin[n + 1) 

Therefore, there is no matter instability. 

For Model pT)) . as it is clear from the identifications ifTO]) . there is no matter instability too. 

In order to study the stability of the de Sitter solution, let us proceed as follows. From the field equations ([2]), we 
obtain the trace 

Df'iR) ^l[R- f'{R)R + V{R) + 't'r] • (27) 

Here, as above, F{R) is F{R) ^ R + f{R) and T = ^'^T^"''. 

Now we consider the (in)stability around the de Sitter solution, where R — Rq, and therefore /(i?o) and f'{Ro), 
are constants. Then since the l.h.s. in Eq. l|27p vanishes for R — Rq, we find 



Ro - /'(i?o)i?o + 2/(i?o) + K^To = . (28) 
Let us expand both sides of l(28|) around i? = i?o as 



R^ Ro + SR. 



(29) 
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One obtains 



Since 



f"{Ro)aSR = ^ (1 - f"{Ro)Ro + f{Ro)) SR . (30) 



(PSR dSR 

UdR = - 

in the de Sitter background, if 



OSR = -^-3Ho—, (31) 



= fi^l 7%R^ > ' ^^2) 

the de Sitter background is stable but, if C(i?o) < 0, the de Sitter background is unstable. The expression for C(i?o) 
could be valid even if f"{Ro) = 0. More precisely, the solution of l(30|) is given by 

SR = A+c^+' + A^c^-* . (33) 

Here A± are constants and 



^ -ago ± y/ggg - cjRo) ^^^^ 

Then, if C{Rq) < 0, A+ is always positive and the perturbation grows up. This leads to the instability. We have also 
to note that, when C{Rq) is positive, if C{Rq) > 9Hq, SR oscillates and the amplitude becomes exponentially small 
being: 

SR= (A cosiUot + B sin ujot)e-''"«'/\ ^ = V^IMElM , (35) 

Here A and B are constant. On the other hand, if C{Ro) < 9Hq, there is no oscillation in SR. 

Let us now consider the case where the matter contribution T can be neglected in the de Sitter background and 
assume f'{R) = in the same background. We can assume that there are two de Sitter background solutions satisfying 
/'(i?) — 0, for R — Ri and R — R2 as it could be the physical case if one asks for an inflationary and a dark energy 
epoch. We also assume f'{R) 7^ if i?i < R < R2 or R2 < R < Ri. In the case C{Ri) < and C(i?2) > 0, 
the de Sitter solution, corresponding to i? = is unstable but the solution corresponding to R = R2 is stable. 
Then there should be a solution where the (nearly) de Sitter solution corresponding to i?i transits to the (nearly) de 
Sitter solution R2. Since the solution corresponding to R2 is stable, the universe remains in the de Sitter solution 
corresponding to R2 and there is no more transition to any other de Sitter solution. 

As an example, we consider Model I. For large curvature values, we find 

MR) = -^+J^- (36) 
Here A and a are positive constants and n is a positive integer. Then we find 

1 n2n+2 

^(^)-m-M2;rTT)^>«- 

This means that the de Sitter solution in Model I can be stable. We have also to note that C(i?o) ^ Hq"''^* /m^"'^^ . 
Here is the mass scale introduced in [16] and <C HqI this means that C(i?o) ^ QHq and therefore there could 
be no oscillation. 

We may also consider the model proposed in [l^(here Model V): 

fv(R) = —. • 38 

Here a, (3, and 7 are positive constants and n is a positive integer. In Figl2l we show the behavior of Model V and of 
its first derivative. When the curvature is large {R 00), f{R) behaves as a power law. Since the derivative of f{R) 
is given by 

, , , (a7i?2" - 2ai?" - /3) 
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we find that the curvature Rq in the present universe, which satisfies the condition f'{Ro) = 0, is given by 



1 + A 1 



/37 



l/n 



and 



T \ 2+ y/1 + pj/a 

As shown in [l^ , the magnitudes of the parameters is given by 

a - 2i?oi?(7^", P ^4RIR^^"R'^-\ 7-2J?o-R(7^"i?r^- 
Here Rj is the curvature in the inflationary epoch and we have assumed f{Ri) ^ {a/'^)Rf ^ Rj. 



(40) 



(41) 



(42) 



£ 0.5 




Figure 2: Plots of Model V l[38|l (solid line) and its first derivative (daslied line). Here n = 2 and a,/?, 7 are assumed as in 
l|42p witli tlie value of Ro taken in tlie Solar System. f'{R) is negative for < 7i < 0.64. f{R) is given in tlie range < iZ < 1 
where we have adopted suitable units. 



C(i?o) in dMl) is given by 



By using the relations (|42|. we flnd 



C{Ro) 



1 + 7i?S' 



/"(i?o) 2n2ai?2n-2 (7i?» - 1) 



C{Ro)- 



An'^Ro 



(43) 



(44) 



which is positive and therefore the de Sitter solution is stable. We notice that C{Rq) < 9Hq and therefore, there 
could occur oscillations as in l(35|) . 

Furthermore, we can take into account the following model [20] (Model VI): 



fviiR) = -a 



tanh(M^)+tanh(^ 



^b{R-Ro) _ I ^bRo _ I 
Qb{R-Ro) _|_ I QbRo _|. I 



(45) 
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where a and h are positive constants. When i? — > 0, we find that 

ahR 



fvi{R) ^ - 

and thus /(O) = 0. On the other hand, when R +oo, 



2cosh2 (^) 



fvi{R) -2Aeff = -a 



1 + tanh 



fbRo 



If Rq, in the present universe, Aofr plays the role of the effective cosmological constant. We also obtain 

ab 



fvAR) = 



which has a minimum when R = Rn, that is: 



2cosh2(^(^) 
fvi{Ro) = - 



ab 



Then in order to avoid anti-gravity, we find 



o<i + /(,,(i?o) = i- Y 



(46) 



(47) 



(48) 



(49) 



(50) 



Beside the above model, we can consider a model which is able to describe, in principle, both the early infiation 
and the late acceleration epochs. The following two-step model (Model VII): 



fvii{R) = -"0 



tanh 



&o {R - Ro) 



tanh 



boRo 



tanh 



bi (i? - Ri) 



tanh 



biRi 



could be useful to this goal. Let us assume 

Ri :S> Ro , a/ » ao , 6/ < bo 

and 

biRi > 1 . 

When i? ^ or i? < i?o < -R/, fvii{R) behaves as 

ao&o , aibi 



fvii{R) 



2cosh2(^) 2cosh2(^) 
and we find again /y//(0) — 0. When R^ Rj, we find 



R 



fiR) 



VII 



-2A, 



-ao 



1 + tanh 



boRo 



1 + tanh 



biRi 



1 + tanh 



biRi 



On the other hand, when Ro ^ R^ Ri, we find 

6oi?o 



fvii{R) — > -ao 



1 + tanh 



ajbjR 



2cosh2 (^) 



-2Ao = —ao 1 + tanh 



boRo 



Here, we have assumed the condition l(53|) . We also find 



fviiiR) 



ajbj 



2 ( bo(R-Ro) 



2 cosh 



2 cosh^ ( 2 

which has two minima for R ^ Ro and R ^ Ri. When R — Rq, we obtain 

a lb I 



1 ( bi{R~B.i) \ ' 



fviiiRo) = -"0^0 



2cosh2(M«i^) 



> -ajbj - aobo . 



(51) 

(52) 
(53) 

(54) 

(55) 
(56) 
(57) 

(58) 
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On the other hand, when R = Rj, we get 

fviARi) = -aibi 



ao^o 



> —aibi — aobo 



Then, in order to avoid the anti-gravity behavior, we find 

ajbi + aobo < 1 



(59) 



(60) 



Let us now investigate the correction to the Newton potential and the matter instabiUty issue related to Models VI 
and VII. In the Solar System domain, on or inside the Earth, where R » i?o, f{R) in Eq.l(45|) can be approximated 

by 

fvi{R} ^ -2Aeff + 2ae-^(«-««) . (61) 
On the other hand, since i?o <C <C i?/ , by assuming Eq. ijSS]). f{R) in can be also approximated by 

fviiiR) - -2Ao + 2ae-'"'(^-«o) ^ (62) 



which has the same expression, after having identified Aq = AcS and bo — b. Then, we may check the case of 
only. In this case, the effective mass has the following form 



jb(i?-i?,o) 



(63) 



which could be again very large. In fact, in the Solar System, we find R - IQ-^^eY^. Even if we choose a ~ 
1/6 ^ Rq ^ (lO^'^^eV) , we find that ^ 10^''"'°eV^, which is, ultimately, extremely heavy. Then, there will 
be no appreciable correction to the Newton law. In the Earth atmosphere, R ^ 10^^" eV^, and even if we choose 



10,000,000,000 ^\r2 



eV . Then, a correction to the Newton law is 



a ~ 1/& ~ i?o ~ (10 ^^eV) again, we find that ml ~ 10 
never observed in such models. In this case, we find that the effective potential U{Rb) has the form 

t/(i?e) = -:^('2A+iV-^^"^ 



2ab 

which could be negative, what would suppress any instability. 

In order that a de Sitter solution exists in /(i?) -gravity, the following condition has to be satisfied: 

R - RfiR) - 2/(i?) . 

For the model l(45|) . the r.h.s of l(65|) has the following form: 



R=-- 



baR 



2coshMMii_iM 



2a 



tanh 



biR-Ro) 



tanh 



bRo 
2 



(64) 



(65) 



(66) 



For large R, the r.h.s. behaves as 

baR 



2 cosh' 



2a 



tanh 



b {R - Ro) 



tanh 



bRo 
2 



2a 



(67) 



although the l.h.s. goes to infinity. On the other hand, when R is small, the r.h.s. behaves as 



baR 



2 cosh' 



+ 2a 



tanh 



b{R- Ro 



tanh 



bRo 
2 



baR 



2cosh2(Mo) 



Then if 



ba 

2cosh^ (^) 



> 1 



(68) 



(69) 
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there is a de Sitter solution. Combining Eg . ((691) with Eq. ljSOp . we find 

2 > ab> ^ 



2cosh^ (^) 



The stability, as above, is given by C(i?ds), where Rds is the solution of (|66| . The expression is given by 

2cosh3(^fc^^ 



C{R. 



-Rds 



ab"^ sinh ^ 
Let us now rewrite Eq. ljMl) as follows, 

'6(i?ds-i?o) 



dS 



2a 



tanh 



tanh 



2 



b tanh 



1 + 



1 

b(-RdS--Ro) 
2 



(70) 



(71) 



ab 



2 cosh 



2 ^ b(fldS--Ro) ^ 



(72) 



Then by using l(72|). we may rewrite iffTI) in the following form: 

~a%^ (1 - x2) \{x - xof + 



where 



and therefore we have 



C{Rds) 



X = tanh 



ab'^x (1 - x^) [2 + ab{l- x^ 
b {Rds — Ro 



Xq 



tanh 



bRo 
2 



(73) 



(74) 



1 < Xq < X < 1 , Xq < 



(75) 

Let us now consider ((66l) in order to find a de Sitter solution. Since Eq. l|66p is difficult to solve in general, we assume 
< Rds < ^0- Then we find 

2cosh2 



RdS — 7 

bxQ 



6=1- 



ab 



= 1 - 



ab (1 — a;o) ' 



(76) 



Eq. l(69l) tells that the parameter e is positive and, by assumption, very small: < e <C 1. Since e is small, by using 
Eqs.jlll), we find 



X = Xq + 



2x0 



Then by using the expression (|73|) for C{Rds), we find 

C(i?ds' 

From the definition of e in l(76|) , we find 



ab^xo (1 - xl) [2 + ab{l~ xl)] 



(77) 

(78) 

(79) 
(80) 



ab{l-xl) =2 + 2e + 0{e^) , 
and then, from Eq.lfTO]). Eq. ((78|) can be written as follows; 

C(i?ds)~-T^. 

bxo 

Since a:o < in the condition l(75l) . we find C{Rds) > and therefore the de Sitter solution is stable. 

In Fig. [3l we have plotted the two models l(45|) and fSTj) written in the form F{R) = R + f{R)- We have 
used the inequalities (f52|) assuming, Rj ^ pg ^ 10^^'*g/cm'^ for the Galactic density in the Solar vicinity and 
i?o ~ Pg ~ 10~^^g/cm'^ for the present cosmological density. . 

Our task is now to find reliable experimental bounds for such models working at small and large scales. To this 
goal, we shall take into account constraints coming from Solar System experiments (which, at present, are capable of 
giving upper limits on the PPN parameters) and constraints coming from interferometers, in particular those giving 
limits on the (eventual) scalar components of GWs. If constraints (and in particular the ranges of model parameters 
given by them) are comparable, this could constitute, besides other experimental and observational probes, a good 
hint to achieve a self-consistent f{R) theory at very different scales. 
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Figure 3: Plots of Model VI (gS]) (solid line) and Model VII (HU (dashed line). Here 6 = 2 and 6/ = 0.5 with a = 1.5 and 
ai — 2. The value of Ri is taken in the Solar System while Ro corresponds to the present cosmological value. 



Mercury perihelion Shift 


|27-/3- 1| < 3 X 10~^ 


Lunar Laser Ranging 


4/3-7-3 = (0.7 ±1) X 10-^ 


Very Long Baseline Interferometer 


|7 - 1| < 4 X 10"* 


Cassini Spacecraft 


7 - 1 = (2.1 ± 2.3) X 10"^ 



Table I: Solar System experimental constraints on the PPN parameters. 



IV. CONSTRAINING /(i?)-MODELS BY PPN PARAMETERS 

The above models can be constrained at Solar System level by considering the PPN formaHsm. This approach is 
extremely important in order to test gravitational theories and to compare them with GR. As it is shown in [33, 3S], 
one can derive the PPN-parameters 7 and /? in terms of a generic analytic function F{R) and its derivative 

- ^ - ^ ^ (81) 



F'{R) + 2F"{Ry 



f3~l 



F'{R)-F"{R) 



2F'{R)+?,F"{RY 



c?7 



(82) 



These quantities have to fulfill the constraints coming from the Solar Systeni experimental tests summarized in Table 
I. They are the periheHon shift of Mercury [4§|, the Lunar Laser Ranging fs^], the upper Hmits coming from the Very 
Long Baseline Interferometry (VLBI) [5l| and the results obtained from the Cassini spacecraft mission in the delay 
of the radio waves transmission near the Solar conjunction [H^l- 



Let us take into account before the /(i?)-models lfT0l) - lfT3|) . Specifically, we want to investigate the values or the 
ranges of parameters in which they match the Solar-System experimental constraints in Table HI In other words, we 
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use these models to search under what circumstances it is possible to significantly address cosmological observations 
by /(i?)-gravity and, simultaneously, evade the local tests of gravity. 

By integrating Eqs. l(8T|) - (|82|) . one obtains f{R) solutions depending on /3 and 7 which has to be confronted with 
Pexp and jexp [38}. If we plug into such equations the models (fT0l) - lfT3|) and the experimental values of PPN parameters, 
we will obtain algebraic constraints for the phenomenological parameters {n,p,q, X, s}. This is the issue which we 
want to take into account in this section. 

From Eq.dHI]), assuming F'{R) + 2F"(Rf ^ and defining A = 

[F"{R)f - AF'{R) = Q. (83) 

The general solution of such an equation is a polynomial function |38| . 
Considering Model II given by (fTTI) . '^^ obtain 



1-7 
27-1 



we obtain 



2pi? § 



-p-i 



Rc 



7-1 



27-1 



2 / R.^ 



RURl-{2p+l)R^fx^ 



[R^ + RlT 



0. 



(84) 



Our issue is now to find the values of A, p, and R/Rc for which the Solar System experimental constraints are satisfied. 
Some preliminary considerations are in order at this point. Considering the de Sitter solution achieved from (fTT]) . we 
have R — const — Ri — xiRc, and xi > 0. It is straightforward to obtain 



xi (l + xj) 



2\P+1 



'i+xir+'^i-{p+i)xi 



(85) 



On the other hand, the stability conditions Fn > and Fub > give the inequality 



(1 + xlY^^ >l + (p + 2)xl + {p+l) {2p +l)xt, (86) 

g 

which has to be satisfied. In particular, for p = 1, it is xi > \/3 and then A > — -;= = 1.5396. In addition, the value 

3v 3 

of xi satisfying the relation ((86|) is also the point where A(a;i), in Eq. l|85p . reaches its minimum. 

To determine values of R compatible with PPN constraints, let us consider the trace of the field equations ^ and 
expHcit solutions, given the density profile p{r), in the Solar vicinity. One can set the boundary condition considering 



Fr^^ = F^RiR = epg) (87) 

where pg ^ 10~^''g/cm^ is the observed Galactic density in the Solar neighborhoods. At this point, we can see 
when the relation l(84|) satisfies the constraints for very Long Baseline Interferometer (7 — 1 = 4 x 10""*) and Cassini 
Spacecraft (7 — 1 = 2.1 x 10^5). This allows to find out suitable values for p. 

An important remark is in order at this point. These constraint equations work if stability conditions hold. In the 
range 

0<^< ' 



Rc V2^+T 

F,RR is negative for the model ifTTj) and then stability conditions are violated. To avoid this range, we need, at least, 
> 1. For example, we can choose = 3.38, corresponding to de Sitter behavior. Then we have p = 1 and A = 2. 

On the other hand, for 0.944 < A < 0.966, we have p = 2 and = -s/S; finally, for R >> Rc, we have A = 2 and 
p = 1.5. For these values of parameters, the Solar System tests are evaded. 

Let us consider now Model I, given by Inserting it into the relation l(83|) . we get 
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f'(R) of I model 

f'(R) of II model 

f'(R) of III model 

f'(R) of IV model 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



R/R =x 

c 



Figure 4: Plots of the first derivatives of four different models as function of a; = Model I (dashed) is drawn for n = 1 and 
A = 2. Model II (dashdot), for p = 2, A = 0.95. Model III (dotted), for s = 0.5 and A = 1.5. Model IV (solid) is for q = 0.5 
and A = 0.5. The labelled values of x indicate where the derivative changes its sign. 



i?3 



(JL 



2n 



R 



2n 



/ \ 271 

2n(f) i?,A 


7-1 
27-1 




(27i + l)(^-£j -2n+l 


i?4 




2n 
) +1 


6 



1 2 



N in 
R \ td2 \2 
"R7 ) 



= 



Using the same procedure as above, A is related to the de Sitter behavior. This means 



(89) 



A 



(1 + xf 



a;2"-i (2 + 2xf - 2n) ' 
while, from the stability conditions, we get 

2xt - (271 - 1) (2n + 4) x?" + (2?! - 1) (2n - 2) > . 
For n — 1, one obtains xi > -\/3 , A > In this model, F_Rfl is negative for 

^ R f2n-l\^ 

The VLBI constraint is satisfied for n = 1 and A = 2, while, for n = 1 and A = 1.5, Cassini constraint holds. 
By inserting Model III, given by Eq. lfT2l) . into the relation (|83l) . we obtain 



(90) 



(91) 



(92) 



i?-2si?e(^)"A 





7-1 




27-1 



R^ 



(93) 
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Figure 5: As above for the second derivatives of the models. 



The de-Sitter point corresponds to 

^23 + 1 



2(a;f - s - 1) 



(94) 



while the stabihty condition is > 2s^ + 3s + 1. VLBI and Cassini constraints are satisfied by the sets of values: 



s = 1, A = 1.53, for ~ 1; s = 2, A = 0.95, for = \/3, ; s = 1, A = 2, for = 3.38. 



Finally let us consider Model VI, given by Eq.lgl]), and Model VII, given by Eq.lHTj). Using Eq.(l83l) for (jUl), we 

get 



^6asech^ (^b{R — Rq 



b-'asech^ -b{R - Rq) tanh^ -b{R - Rq) 



1. 



7-1 



27- 1 



(95) 



As above, considering the stability conditions and the de Sitter behavior, we get the parameter ranges < 6 < 2 and 
< a < 2 which satisfy both VLBI and Cassini constraints. Inserting now Model VII in l(83|l . we have 



7-1 



27-1 
1 



6asech^ ( h{R - i?o) j - 6/a/sech^ (^^^^^ -Ri)] +2 



b^asech^ (^^b{R - Ro)^ tanh - ^0)^ - bjaisech^ {^^^'^^ ~ ^^)) {\^'^^ ^ 



2 (96) 



From the stability condition, we have that F^n > for R > 0, (see Figl6]) and F^nn < for < i? < 2.35 in suitable 
units (see Figl?]). Observational constraints from VLBI and Cassini experiments are fulfilled for 



(97) 



Plots for b = 2, bj = 0.5, a — 1.5 and aj = 2, verifying the constraints, are reported in Figs. [6] and El 
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2 - . 



-f'(R) of VI model 
-f'(R) of VII model 




Figure 6: Plots represent the first derivatives of functions l|50p (solid line) and l|5ip (dashed line). Here, 6 = 2, 6/ = 0.5, a = 1.5 
and ai — 2 with Ri with the Solar System value and Ro the today cosmological value. It is F^R > for i? > 0. 



Considering now the relation for (3 given by Eq. i(82|) . one can easily verify that it is 



dR 



d_ 

dR 



F"{Rf 



F'{R) + 2F"{RY 



= 0, 



and this result implies 



4(/3-l) = 0. 



(98) 



(99) 



This means the complete compatibility of the f{R) solutions between the PPN-parameters (3 and 7. 

Now we want to see if the parameter values, obtained for these models, are compatible with bounds coming from 
the stochastic background of GWs achieved by interferometric experiments. 

V. STOCHASTIC BACKGROUNDS OF GRAVITATIONAL WAVES TO CONSTRAIN /(7?)-GRAVITY 



As we said before, also the stochastic background of GWs can be taken into account in order to constrain models. 
This approach could reveal very interesting because production of primordial GWs could be a robust prediction for 
any model attempting to describe the cosmological evolution at primordial epochs. However, bursts of gravitational 
radiation emitted from a large number of unresolved and uncorrelated astrophysical sources generate a stochastic 
background at more recent epochs, immediately following the onset of galaxy formation. Thus, astrophysical back- 
grounds might overwhelm the primordial one and their investigation provides important constraints on the signal 
detectability coming from the very early Universe, up to the bounds of the Planck epoch and the initial singularity 
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Figure 7: Second derivatives of Model VI (solid line) and VII (dashed line). Here F^rr is negative in the range < J? < 4 for 
Model VI and in the range Q < R < 2.35 for Model VII. As above, we have used b = 2, bi = 0.5, a = 1.5 and ai = 2 with the 
value of Ri taken in the Solar System and Ro for the today cosmological value. 



It is worth stressing the unavoidable and fundamental character of such a mechanism. It directly derives from the 
inflationary scenario [5l,[5§|, which well fits the WMAP data with particular good agreement with almost exponential 
inflation and spectral index « 1, [sol. Isij. 

The main characteristics of the gravitational backgrounds produced by cosmological sources depend both on the 
emission properties of each single source and on the source rate evolution with redshift. It is therefore interesting to 
compare and contrast the probing power of these classes of /(i?)-models at hight, intermediate and zero redshift [6a |. 

To this purpose, let us take into account the primordial physical process which gave rise to a characteristic spectrum 
i^sgw for the early stochastic background of relic scalar GWs by which we can recast the further degrees of freedom 
coming from fourth-order gravity. This approach can greatly contribute to constrain viable cosmological models. The 
physical process related to the production has been analyzed, for example, in [H, HI, HH] but only for the first two 
tensorial components due to standard General Relativity. Actually the process can be improved considering also the 
third scalar-tensor component strictly related to the further f{R) degrees of freedom [i^ . 

Before starting with the analysis, it has to be emphasized that the stochastic background of scalar GWs can be 
described in terms of a scalar field $ and characterized by a dimensionless spectrum (see the analogous definitions 
for tensorial waves in ^45.. ,54. .Sa. ,57,]). We can write the energy density of scalar GWs in terms of the closure energy 
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density of GWs per logarithmic frequency interval as 

ns,Uf)--^, (100) 

Pc dlnf 

where 



^ 8^ 

is the critical energy density of the Universe, Hq the today observed Hubble expansion rate, and dpsgw is the energy 
density of the gravitational radiation scalar part contained in the frequency range from f to f + df. We are considering 
now standard units. 

The calculation for a simple inflationary model can be performed assuming that the early Universe is described by 
an inflationary de Sitter phase emerging in the radiation dominated era [El. Issl. [57| . The conformal metric element is 

ds^ = a^{Tj)[-dri^ + d'x^ + hf,^{r], ~x )dx^'dx''], (102) 

and a GW with tensor and scalar modes in the z+ direction is given by [47] 

V(i - ^) = A+{t - z)eW +A-{t- z)e|,^,) + $(t - z)e(fj . (103) 

The pure scalar component is then 

V = ^e^^il , (104) 

where e^^J is the polarization tensor. 

It is possible to write an expression for the energy density of the stochastic relic scalar gravitons in the frequency 
interval (w,a; + duj) as 

where /, as above, is the frequency in the standard comoving time. Eq. l|105p can be written in terms of the today 
and de Sitter values of energy density being 

Ho = , Hds - o„2 ■ (106) 



Introducing the Planck density ppianck — "^t^, the spectrum is given by 



3c2 ' 3c2 



^ _ 1 dpsgn, _ f dpsgw _ 16 Pds I . 

Pc din J Pc df 9 ppianck 

At this point, some comments are in order. First of all, such a calculation works for a simplified model which does 
not include the matter dominated era. If such an era is also included, the redshift at equivalence epoch has to be 
considered. Taking into account also results in [HHl, we get 

^sg^{f) = ^^^{l + Zc,)-\ (108) 
y PPlanck 

for the waves which, at the epoch in which the Universe becomes matter dominated, have a frequency higher than 
Hcq, the Hubble parameter at equivalence. This situation corresponds to frequencies / > (1 + Zerj)^^^-ffo- The redshift 
correction in Eq. l|108p is needed since the today observed Hubble parameter Hn would result different without a 
matter dominated contribution. At lower frequencies, the spectrum is given by [m . IsH I 

nsgn,{f) . (109) 

Nevertheless, since the spectrum falls oflf cx /^^ at low frequencies, this means that today, at LIGO- VIRGO and LISA 
frequencies (indicated in Fig. [8]), one gets 

^^s3»(/)^?oo < 2.3x 10-1^ (110) 
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^dWOHz) < 2 X lO"^'^ 


LIGO 


^dWOHz) < 2 X 10-2^ 


VIRGO 


^c{WOHz) < 2 X 10"^' 


LISA 



Table II: Upper limits on the expected amplitude for the GW scalar component. 



It is interesting to calculate the corresponding strain at « lOOHz, where interferometers like VIRGO and LIGO reach 
a maximum in sensitivity (see e.g. |63 . [63|). The well known equation for the characteristic amplitude [H, [55| . 
adapted to the scalar component of GWs, can be used. It is 

^ 1.26 X 10-i« (^i^) ^hl,,n,gUf ), (111) 
and then we obtain the values in the Table Ull 



Energy 



LISA LIGO 




Hz ( - Graphics 



Figure 8: The spectrum of relic scalar GWs in inflationary models is flat over a wide range of frequencies. The horizontal axis 
is logjQ of frequency, in Hz. The vertical axis is logj^Q ilgsw The inflationary spectrum rises quickly at low frequencies (wave 
which re-entered in the Hubble sphere after the Universe became matter dominated) and falls off above the (appropriately 
redshifted) frequency scale fmax associated with the fastest characteristic time of the phase transition at the end of inflation. 
The amplitude of the flat region depends only on the energy density during the inflationary stage; we have chosen the largest 
amplitude consistent with the WMAP constrains on scalar perturbations. This means that, at LIGO and LISA frequencies, we 
have ttsg^ < 2.3 * 10"^^ 



In summary, the above results point out that a further scalar component of GWs, coming e.g. from /(i?)-gravity, 
should be seriously considered in the signal detection of interferometers. As discussed in [63|, this fact could constitute 
either an independent test for alternative theories of gravity or a further probe of GR capable of ruling out other 
theories. 

At this point, using the above LIGO, VIRGO and LISA upper bounds, calculated for the characteristic amplitude 
of GW scalar component, let us test the /(i?)-gravity models, considered in the previous sections, to see whether they 
are compatible both with the Solar System and GW stochastic background. 
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Before starting with the analysis, taking into account the discussion in Sec. II, we have that the GW scalar component 
is derived considering 



$ = -- 



5o_ 



a = -Hl + f{A)) = \nF'{A), 5a=J^^^5A. 



As standard, we are assuming small perturbations in the conformal frame [47|)- This means 



(112) 



(113) 



These assumptions allow to derive the "linearized" curvature invariants R^vpa , Rfj.i' and R and then the linearized 
field equations [Hi 



(114) 



As above, for the considered models, we have to determine the values of the characteristic parameters which are 
compatible with both Solar System and GW stochastic background. 

Let us start, for example, with the model (fT2|) . Starting from the definitions l|112p . it is straightforward to derive 
the scalar component amplitude 



^iii 



,s(2s + l)(f)^^+^ 



si?c(f) A-i? log 2-2s(i|) 



(115) 



Such an equation satisfies the constraints in Table ITTl for the values s — 0.5, '--^ 1, A = 1.53 and s = 1, ~ 1, 

A 0.95 (LIGO); s = 2, ;^ = V3, A==2 (VIRGO); s = 1, A = 2 and = 3.38 (LISA). 

It is important to stress the nice agreement with the figures achieved from the PPN constraints. In this case, we 
have assumed i?c ~ Pc ~ 10~2^g/cm'^, where pc is the present day cosmological density. 

Considering the model ([9]), we obtain 



n 


(2n+l)(^-£j -2n+l 




^ 2n-l 

) ' 


(f) +1 


\ 




(f) +1 


i?,A 





(116) 



The expected constraints for GW scalar amplitude are fulfilled for n — 1 and A = 2 and for n = 1 and A = 1.5 when 
0.3 < ;^ < 1. 

Furthermore, considering the model l|lip . one gets 



2p (l + 1^) " Rc ((1 + 2p) i?2 - Rl) X 






In 









(117) 



The LIGO upper bound is fulfilled for p = 1, > V3, A > the VIRGO one for p = 1, = 3.38, A = 2; finally, 

for LISA, we have p — 2, ^ — \/3 and 0.944 < A < 0.966. Besides, considering LISA in the regime R >> Rc, we 
have A = 2 and p — 1.5. 

Finally, let us consider Models VI and VII. We have 



VI 



b^at&nh [^b{R- Ro)] 



[ba + cosh{b{R - Rq)) + 1] In 



ba 



cosh(b(_R-i?.o)) + l 



(118) 
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and 

<^yjj = log [0.5 (6Q;sech^(0.56(i? - i?o)) - 5/a/sech^(0.5fe/(i? - i?/)) + 2)] 

X [6asech2(0.56(i? - i?o)) - biaisech^iOMiiR - i?/)) + 4] (119) 
X [6^asech^(0.56(i? - i?o)) tanh(0.56(i? - Ro)) - bj a i sech^{0. 5b i{R - i?/)) tanh(0.56/(i? - i?/))] . 

These equations satisfy the constraints for VIRGO, LIGO and LISA for b = 2, bj = 0.5, a = 1.5 and aj = 2 with 
Ri valued at Solar System scale and Rq at cosmological scale. 

VI. CONCLUSIONS 

In this paper, we have investigated the possibility that some viable f{R) models could be constrained considering 
both Solar System experiments and upper bounds on the stochastic background of gravitational radiation. Such 
bounds come from interferometric ground-based (VIRGO and LIGO) and space (LISA) experiments. The underlying 
philosophy is to show that the f{R) approach, in order to describe consistently the observed universe, should be tested 
at very different scales, that is at very different redshifts. In other words, such a proposal could partially contribute 
to remove the unpleasant degeneracy affecting the wide class of dark energy models, today on the ground. 

Beside the request to evade the Solar System tests, new methods have been recently proposed to investigate the 
evolution and the power spectrum of cosmological perturbations in f{R) models [131 ■ The investigation of stochastic 
background, in particular of the scalar component of GWs coming from the f{R) additional degrees of freedom, could 
acquire, if revealed by the running and forthcoming experiments, a fundamental importance to discriminate among 
the various gravity theories (gII- These data (today only upper bounds coming from simulations) if combined with 
Solar System tests, CMBR anisotropics, LSS, etc. could greatly help to achieve a self-consistent cosmology bypassing 
the shortcomings of ACDM model. 

Specifically, we have taken into account some broken power law f{R) models fulfilling the main cosmological 
requirements which are to match the today observed accelerated expansion and the correct behavior in early epochs. 
In principle, the adopted parameterization allows to fit data at extragalactic and cosmological scales [l^. Furthermore, 
such models are constructed to evade the Solar System experimental tests. Beside these broken power laws, we have 
considered also two models capable of reproducing the effective cosmological constant, the early infiation and the late 
acceleration epochs jl^l- These /(i?) -functions are combinations of hyperbolic tangents. 

We have discussed the behavior of all the considered models. In particular, the problem of stability has been 
addressed determining suitable and physically consistent ranges of parameters. Then we have taken into account 
the results of the main Solar System current experiments. Such results give upper limits on the PPN parameters 
which any self-consistent theory of gravity should satisfy at local scales. Starting from these, we have selected the 
f{R) parameters fulfilling the tests. As a general remark, all the functional forms chosen for f{R) present sets of 
parameters capable of matching the two main PPN quantities, that is jexp and Pexp- This means that, in principle, 
extensions of GR are not a priori excluded as reasonable candidates for gravity theories. To construct such extensions, 
the reconstruction method developed in may be applied. 

The interesting feature, and the main result of this paper, is that such sets of parameters are not in confiict with 
bounds coming from the cosmological stochastic background of GWs. In particular, some sets of parameters reproduce 
quite well both the PPN upper limits and the constraints on the scalar component amplitude of GWs. 

Far to be definitive, these preliminary results indicate that self-consistent models could be achieved comparing 
experimental data at very different scales without extrapolating results obtained only at a given scale. 
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